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Abstract 

Using a holographic proposal for the entanglement entropy we study its be- 
havior in various supergravity backgrounds. We are particularly interested in 
the possibility of using the entanglement entropy as way to detect transitions 
induced by the presence horizons. We consider several geometries with horizons: 
the black hole in AdS^, nonextremal Dp-branes, dyonic black holes asymptot- 
ically to AdSii and also Schwarzschild black holes in global AdSp coordinates. 
Generically, we find that the entanglement entropy does not exhibit a transition, 
that is, one of the two possible configurations always dominates. 



1 Introduction 



Given a system in a pure quantum state |\E' > and density matrix p = >< \[^|, 
if we split the system into two subsystems A and B, the reduced density matrix is 
obtained by tracing over the degrees of freedom in the complementary subsystem, say, 
Pa = TtbP- The entanglement entropy is defined as the von Neumann entropy 

Sa = -TrpAhgpA- (1.1) 

This provides a measure of how entangled or "quantum" a system is. Entanglement 
plays a central role in quantum information theory as it determines the ability to send 
quantum information [1]. The entanglement entropy also plays an important role in 
the study of strongly correlated quantum systems [2]. 

The above definition is completely field theoretical. Interestingly in the context 
where such field theories have supergravity duals, a prescription for the holographic 
computation of the entanglement entropy has been provided [3], that is, a prescription 
for computing Sa completely within the holographic gravity dual. Inspired by the 
Bekenstein-Hawking entropy, it suggests to calculate the entanglement entropy as the 
area associated to a minimal surface '-/a whose boundary is the region A in the field 
theory living at the boundary: 

= (1.2) 

This recipe has been successfully applied to various systems and extended in different 
directions [4] including its covariant formulation [5] . A slightly modified version of the 
entanglement entropy is 

= 7^ / V^. (L3) 

The entropy is obtained by minimizing the above action over all surfaces that approach 
the boundary of the subsystem A. In a very interesting paper [6], it was suggested that 
in the presence of regions with collapsing cycles, which are typical for supergravity 
duals of confining field theories, alternative surfaces arise (see also [7]). By comparing 
the entropy due to two different configurations, it was shown that the entanglement 
entropy could be an order parameter for the confinement /deconfinement transition. 
The motivation comes from the fact that the entanglement entropy jumps from a 
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configuration with dominant term of tlie form A^^ to a configuration witli leading term 

for supergravity backgrounds describing tlieories witli SU{N) gauge groups in tfie 
large and fixed 't Hooft coupling limit. 

One intriguing fact about the results of [6] is that, by analyzing a surface in the 
supergravity dual to the confined phase of the field theory, one is able to anticipate 
the existence of a deconfined phase. This prompts the natural question of whether 
the exploration of the deconfined phase can equally well give information about the 
existence of a confined phase. 

According to the AdS/CFT correspondence [8], the dual of field theories at finite 
temperature involves black hole horizons on the supergravity side [9]. The horizon 
provides a natural end of the space similar to the situation discussed in [6]. In fact, 
in the context of the AdS/CFT is has already been established in various situations 
that there are phase transitions associated to different behavior of surfaces describing 
branes in the presence of a horizon [10-15]. 

In this paper we explore the entanglement entropy for various field theories at 
finite temperature using its holographic definition. Namely, we study the entanglement 
entropy for various black hole geometries. We define our subsystem A to be roughly 
determined by an interval of length / on a curve generated by a spacelike killing vector 
in the conformal boundary of the geometry. Generically, there are two surfaces that 
satisfy the boundary conditions: a smooth one and a piece- wise smooth (see figured]). 
We study the behavior of these two surfaces as a function of the distance /. 

In section |2] we discuss the gravity dual of 2D CFT at finite temperature - the BTZ 
black hole - where the smooth surface is just the geodesic anchored on A. The piece- 
wise smooth surface is the curve composed of the lines stretching from the conformal 
boundary to the horizon connected by the segment of length / on the horizon. We 
show that the smooth configuration is always dominant from the point of view of the 
thermodynamical comparison. Section [3] presents a general setup for the computations 
at hand and discusses explicitly the case of the the supergravity backgrounds describing 
the Af = A plasma and other field theories at finite temperature. Also in section [3] we 
extend the analysis to black p-branes corresponding to various gauge theories on the 
world volume of p-branes at finite temperature. Using the supergravity backgrounds 
for nonextremal Dp-branes for all values of p, we found that the entanglement entropy 
is given by the smooth surface for all /, except for p = 6. In the p = 6 case we observe 
that the entanglement entropy is given by the smooth surface at large I and by the 
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Figure 1: Two competing configurations for the entanglement entropy in tlie presence 
of a black hole horizon. The green surface represents a continuous configuration while 
the red surface goes straight down from the boundary to the horizon. The subsystem 
A is given by the blue section. Its characteristic length is /. 

piece-wise smooth surface at small /. Furthermore in these geometries, we observe that 
the radius of the black hole factors as an overall scaling factor. This scaling implies 
that the entanglement entropy has the same N (number of color) dependence for all the 
temperatures. In section [5l we study the entanglement entropy in black hole geometries 
in global AdS^ and AdS4^. The main motivation for our analysis comes from the fact 
that the Hawking-Page phase transition takes place only in global coordinates [16], 
that is, for the nonextremal Dp-brane geometries there is no Hawking-Page phase 
transition. In the context of the AdS/CFT this is interpreted as a transition for the 
dual field theory on a sphere [9]. Since this feature is absent in the Poincare patch we 
are motivated to study the entanglement entropy in global AdS backgrounds. However, 
we find that in global coordinates, the entanglement entropy is given also by the smooth 
surface. 

Our work suggests that starting from the deconfined phase, changing the lenght of 
the subsystem does not allow for a transition into confinement, as opposed to what was 
observed in [6,7]. According with their results, at zero temperature the length of the 
subsystem seems to play the role of temperature in the confined phase, while in our 
case this lenght corresponds possibly to some other thermodynamical quantity. It is 



3 



probably worth mentioning a recent study of the topological entanglement entropy [24] 
where no change in the dependence on was observed. However there is mounting 
evidence that the entanglement entropy is a bona fide thermodynamical quantity and 
its precise meaning should be illuminated through further work. In section [6] we discuss 
some of the directions suggested by our exploration of the entanglement entropy in the 
context of finite temperature. 



2 2D CFT at finite temperature from the BTZ black hole 

We begin this section with a discussion of the holographic entanglement entropy using 
the gravitational background dual to a 2D CFT at finite temperature. Most of the 
calculation was explicitly done in [4], however, we reproduce it here making emphasis 
on the thermodynamical competition between the two configurations, something that 
was not considered in [4]. The relevant geometry holographically describing the field 
theory at finite temperature is the BTZ black hole [17]: 

ds^ = -(r2 - rl)dt' + -^^dr^ + r^rf^^. (2.I) 

r 

The smooth surface is parametrized by t =constant and r = r{(j)). Subsystem A is 
defined as the region given by < < 2itI/L where L is the circumference of the 
boundary. Thus, / is the characteristic size of A. For the BTZ black hole, the smooth 
surface is just the geodesic in the bulk that connects the two boundary points of A. 
The action of this curve is given by 




Ac= / d<f)Jr^ + - (2.2) 



The equation of motion can be integrated to give 

^ = ^/M(^. (2.3) 
This allows us to relate the length in the direction with the minimum of the curve, 

M = ^l^!jL±l±. (2.4) 

L Tj^ r^- r+ 

The temperature is /3 = RL/r^. The gravitational theory on AdS'^ with radius R is 
dual to a 2D CFT living on its boundary with central charge c = ?>R/2G^^ [18]. 
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The length of the geodesic is then 

Areac = 2R\n — sinh ( —nl] (2.5) 
r_|_ V RL / 

where r^o is a UV cutoff. 

The entanglement entropy given by the continuous configuration is: 

Sc = Arettc = In — sinh (^nl) . (2.6) 
4GJJ^ AG^i' r+ \RL J 

As noted in [3], this result is in agreement with the field theory result for a 2D CFT 
at finite temperature [19]: 

5a = ^og f A sinh ^ V (2.7) 



3 \7ia (3 

where a is a ultraviolet cutoff that can be thought of as a lattice spacing. 

The piece- wise smooth surface is given by the parametrization gluing the 3 surfaces 
described as by = 0, r = r+ and = 21x1/ L. The area for this curve is then 

Too 

f R , 2r+7r/ Too 2r:irl 

Aread = 2 / —==dr + ^— = 2R\n^ + ^—. 2.8 

The first term of the piece-wise smooth configuration consists of two lines that go from 
the boundary to the horizon, their contribution to the entropy is: 

Too 

Sm) = -4t2 / -^^dr = In ^. (2.9) 



»■'+ 

By construction, this entropy is independent of the length / of the interval that defines 
the subsystem A. In intrinsic field theoretic terms we could rewrite it as: 

5d(i) = ^ln^. (2.10) 

O Qj 

The study of entanglement entropies in 1+1 system is well developed. In [20] a formula 
for the entanglement entropy of a system of length / was obtained for conformal field 
theories: (c/3) ln(//a). Some extensions of this result have been discussed recently, 
including systems away from criticality [19,21,22]. Systems close to a phase transition 
(large but finite correlation length) are described by massive quantum field theories 
with mass inversely proportional to the correlation length [19]: 

4 '"(!)■ (^■") 
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The correlation length that appears above is considered to be the inverse of the mass 
m = The correlation length ^ and the inverse temperature (3 can be naturally 
identified in our setup completing the matching of the gravity f l2.10p with the field 
theoretic one. Interpreting the second term in the piece-wise smooth configuration is 
more challenging, we simply note that in field theoretic terms the contribution coming 
directly from the horizon takes the form: 5*^(2) = (c/3)(7r///3). 
Now consider the difference in area: 



AA = Ar-A^ = 2Rln 



sinh 



2r.Trl 



RL 

As expected, the difference is UV finite. Setting it to zero we obtain: 



(2.12) 



sinh(x) = where x = . (2-13) 

RL 

This equation has no solutions and it follows that 

A, < A^. (2.14) 

This implies that the entropy is always given by the smooth surface. As noted in [4] the 
answer matches precisely the field theory calculation [22]. More impressive agreement 
has also been found in the context of disconnected segments in the field theory [23]. 



3 Entanglement entropy for nonextremal Dp branes 

The next geometries that we study are static geometries with a spacelike killing vector 
X^^ that commutes with all other killing vectors. Furthermore we assume that these 
geometries have conformal radius r. Thus there exist a frame where the metric is given 

as, 

ds'^ = -T{r)df + f{r)dr^ + g{r)dx^ + K^dy'dy^ (3.1) 

where x is the affine parameter along and are coordinates of the internal sub- 
manifold. In general, this submanifold can have compact and non-compact directions 
and hij can have non-trivial dependence on the coordinates. However, the geometries 
considered here satisfy: 

hijdy'dy^ = hli,{r)du''du^ + hl^{r, e)de'^de'' (3.2) 
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where the m*'s are non-compact coordinates and 6''^'s are compact. The conformal 
boundary is obtained by taking the large r hmit. The region A that we will consider 
is parametrized by: 

1 I 

t = constant and < x < - 

2 2 

where / is the size of the subsystem A. As discussed above, there are two surfaces that 
satisfy the minimal surface condition. The smooth surface is given by, 

t = constant, and r = r{x) with b.c. r(±//2) = Too, (3-3) 

where r^o is the location of the holographic boundary. The second surface is piece-wise 
smooth and parametrized as: 



X 2 



t = constant, and < 



r = Tq (3-4) 



X 2 



where tq is the lower bound of r. Now we can proceed to study the difference in area 
of these surfaces. 

3.1 Two branches of the holographic entropy 

We start by computing the area of the smooth surface. The induced metric is: 

ds^ = g{r){l + ^r'\x))dx^ + K^dy'dy^ (3.5) 
where prime denotes derivative with respect to x. The volume element is then 



V^= W^7W(l + 44^"(^))^HO/^'(r,^) (3.6) 

V 9{n 

where and h"^ are the determinant of the metric for the non-compact and compact 
sub-manifolds, respectively. Since the non-compact sub manifold will contribute an 
infinite volume, we must work with the area density. It is given by 

A = j dxe^^^^g{r){l + ^r'^{x))h^{r) j d''e^h^{r,e) (3.7) 

where is the dilaton. Before we proceed, we define the following quantities: 

^^W^ = e-^'^^/g{r)h^{r) j d'e^h?{r,e) and /J' = (3.8) 
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The area is: 

A = J dx^/H{r)y/l + (3^r'^{x). (3.9) 

Now we want to find the r(x) that minimizes the area. Since the Lagrangian does not 
exphcitly depend on x, we can use the conserve Hamiltonian to integrate the equation 
of motion. It is then: 

E^L-4L,r')^ . (3.10) 

d{r'y ^1 + /3V'2 

The constant E can be obtained by considering the turning point r^, where r' = 0. We 
thus obtain 

E = ^H{n). (3.11) 

Physically, r^, is the minimum of the surface. It corresponds to x = since the surface 
must be symmetric under x —x. The quantity, r,,, can be used to label different sur- 
faces for different values of /. This relationship is obtained by integrating the equation 
of motion one more time, 

/(r.) = !dx = !% = r -7w7^=^d^- (3-12) 

J J r' Jr, \/H{r) - H{r^) 

Similarly, we can obtain the area as an integral over r, 

A,{n) = J dx^/Wi^^/lVWHx) = j ^H{r) 

Jr, ^Hir)-Hiny 

We also compute the area of the piece- wise smooth surface. The induced line elements 
for different segments are 

ds"^ = f{r)dr'^ + hijdy'dy^ for x = (3.14) 
ds'^ = g{rQ)dx^ + hijdy^dy^ for r = rg. (3.15) 

The area of the piece-wise smooth is then 

Ad = 2 (3{r)^/H{r)dr + l^jH{ro). (3.16) 

In what follows it suffices to identify H and (3 from the geometries; the difference in 
area and / are given by; 
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(3{r)H{r) 



'H{r) - H{n] 
(3{r 



1/2 



H(r) 



H{r) 



1/2 



-.dr. 



'H{r) - H{n) 

It should be noted here that when r^, = tq, the difference in area is negative since 



(3.17) 
(3.18) 



AA(ro) 



P{r)H{r) 



1 _ ^ _ A - ^ 
H(r) [ H(r) 



Similarly if /5 — >^ fast enough as r — Too, then 

AA{n ^ Too) ~ -2 / /3(r)/^ < 0. 



1/2- 



< 0.(3.19) 



(3.20) 



ro 



This statement is almost always true. We will see an exception where (3 (x Thus 
we obtain two important results before we even start to look at specific geometries 



AA(n = r-o) < 0, and Ay4(n 



r„o) < 0. 



(3.21) 



Thus, if the difference in area is monotonous, the smooth surface gives the entanglement 
entropy. Now we proceed to study geometries corresponding to holographic duals of 
the Af = 4 plasma at finite temperature and other field theories leaving in Dp-branes 
world volumes and dyonic black hole. 

3.2 Entanglement entropy in the AT = 4 plasma 

The holographic background dual to the A/" = 4 plasma is a stack of nonextremal 
D3 branes. This background has proved to be an interesting playground for finite 
temperature field theories, in particular, it seems to catch some features displayed by 
experiments at the Relativistic Heavy Ion Collider. The corresponding background 
metric is: 



du^ 



u 



[-hdf + dxidx^) + dQj 



with 



Un 



(3.22) 



(3.23) 
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where the parameter uq completely characterizes the temperature of the background: 
uq = irRT. The coordinate u is related to the holographic coordinate r by r = Ru 
where R is the AdS radius. The quantities H and (3 are: 

1 



H{u) = {R^n^Yu'', (3\u) 



— u, 



4 • 



(3.24) 



The difference in area and / are: 
AA 



2Iul 



I = — 



yea 



y'tdy 



^oJy, Viy'-^)iy'-y'*) 



(3.25) 
(3.26) 



where y = u/uq and / = Q^R^. The y coordinate makes it clear that AA is just scales 
with respect to the temperature. We numerically plot AA as function of Iuq in units 
of 2Iul. 
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(a) 



(b) 



Figure 2: Figure 2(a) shows the behavior of / as a function of for constant Uq. 



Figure 2(b) Shows the difference in area as a function of /. These plots show no phase 



transition the entropy in the nonextremal D3 brane background. 

From figure [2] we see that the difference in entropy does not change sign; instead 
it approach 0. A A is also monotonous with respect to / as observed from [2l / is 
monotonous with respect y*; in fact it logarithmically diverges as y^, — 1 as observed 
from 13.261 Thus by 13.211 A A is always negative. 
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3.3 Entanglement entropy for gauge theories on the world volumes of Dp 
branes 



The backgrounds dual to finite temperature field theories are given in the string frame 
as [25]: 



u 



1- 



u 



-dx^ 



+ 



-1 



u 



a/2 



u 



du^+^u^'^dnl 



8-p 



where a = 7 — p and b = p — 3. The dilaton is also given as: 

-b/2 



u 



a/2 



(3.27) 



(3.28) 



It is important to note that for p < 3, the radius of the (8 — p)-sphere diverges in the 
UV limit and shrinks in the IR limit. However, for p > 3, it shrinks to zero in the 
UV limit and expands in the IR limit. At p = 3, it decouples from the AdS sector. A 
similar feature takes place for the dilaton. The quantities H and (3 are 

^2 



A 



where 



C 



(27r) 



2(p-2) 



9ym 



A 



2(9-p)- 



(3.29) 



(3.30) 



The difference in areas AA and / as integrals over y are 
AA = A,- Ad 



2Cul 



I 



.a-2 



y2a_ya^y2a_y2a_yc 

y'^dy 



■dy 



y" dy 



^^y2a_y2a) (ya _ 1) 



(3.31) 



We observe something interesting here. The difference in areas AA scales as Uq for all 

p — b 

p while / scales as . For p > 5, it is proportional to uq. In figure [3] we plot the 

P-5 

difference in areas in units of Cuq with respect to / in units of v Amq^ for p = 3, 4, 5. 
In figure H] we show the plot for p = 6. 

Here we observe a transition for p = 6 and no transition for p < 6. This result is 
also in agreement with result 13.211 
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Figure 3: Figure 3(a) shows the behavior of /M^f ^^^"^ as a function of at p = 3,4,5 
(Blue, Green, Orange respectively). Figure 3(b) Shows the difference in area as a 



function of /. These plots show no phase transition for the entanglment entropy in the 
backgrounds of nonextremal D3, D4 and D5 branes. 
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Figure 4: Figure |l](a) shows the behavior of as a function of for constant uq 

at p = 6. Figure 111(b) Shows the difference in area as a function of /. These plots show 
a transition for the entanglement entropy in the nonextremal D6 brackground. 



4 Dyonic Black Hole 

This is a very important system from the holographic point of view. It has the potential 
to describe a few systems that are certainly of interest from the condensed matter point 
of view. A partial list of interesting applications can be found in [28]. 
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4.1 The solution 

The dyonic black hole is a solution to Einstein Maxwell on AdS^. The solution is a 
consistent truncation of 11 dimensional supergravity on AdS^ x S'^. Some interesting 
properties of this solution and its potential applications to condensed matter have been 
recently discussed in various papers including, [26-28] 
The relevant metric is: 

p2 J 2 

ds' = R^u\-h{u)dt' + dz" + dy^) + (4.1) 



h{u) 



with 



4 3 



h{u) = l + {h' + g^)^ - (1 + /.^ + g^)^, (4.2) 
and electromagnetic field tensor 

F = huldx Ady - quldt A du~^ . (4.3) 

The three parameters mq, h and q are related to the physical quantities, mass, electric 
charge and magnetic field of the black hole in the following way, 

M = a(l + /i^ + q^)ul, C = qui, and, B = hul (4.4) 

respectively. The quantity a is a constant independent of h, q and Uq. It is useful to 
also define the parameters = h"^ + q^ and = p^u^ to characterize the effect of the 
electromagnetic field. In terms of the physical parameters, uq and are given by 

n2 

= ^, Mo + - "^^0 = (4.5) 

where m = M/a. A computation of the mass using holographic renormalization can 
be found in [27]. Equation 14.51 has a positive real root for uq when 



4 



which implies 



< ( ^ ) • (4-7) 



From these relationship, we can expect extremality at p^ = 3. Furthermore, we observe 
that even though Q is bounded by the mass, p^ is allowed to take on any values. Thus 
if there is extremality at = 3, we can expect different thermodynamical description 
of space when < 3 from when > 3. This is discussed in the next section. 
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4.2 Comments on thermodynamics 

With the coordinate change y = u/uq, h{u) can be rewritten 

Ky) = -Ay' + p'-ii + p')y) = -Ay'-y + p'ii-y)) 
y y 

= '^{y' + y' + y-p')^^J(y) 

= (,^H-(l+,o)i/+^) (4.8) 

where po satisfies 

pI + pI + Po-p' = 0. (4.9) 
We observe that h{y) has at least one positive zero at y = 1 corresponding to u = uq; 

and at most two positive zeros with y = po corresponding to m = PqUq. We can 
understand this second zero by studying equation 04.91) . We can rewrite it as 

p'(Po) = Po + Po + Po. (4.10) 

First we observe that po vanishes when p^ = and po = 1 when p^ = 3. In addition, 
the function p^(po) is monotonic since its derivative has no real zeros. Thus p^ has a 
one to one relationship with po. This justifies the above statement that h{u) has at 
most two real zeros. Since p^ is not bounded, po is also unbounded. The inverse of this 
function po(p^) is then: 

1/3 

■ (4.11) 



7 + 27p^ + 3V3 V3 + 14p2 + 27p^ 



¥-h-2 , 

Po = ^7 where h -- 

36 

The following picture emerges for the position of the horizon p*Mo. 

{1 when < p2 < 3 
(4.12) 
Po when 3 < p < oo 

The temperature of the horizon is given by the surface gravity. 



K 1 J-(V'(f)(S7„i^)\,„„„ 



^-T.-T.i (4.13) 



where is a timelike killing vector. With = (^)", the temperature is. 



no 



3 - p2 for < p2 < 3 



T = ^{ - . • (4.14) 



47r 1 2pg+po+p"(po-2) 



3 for 3 < p < oo 



We observe extremality at p^ = 3 as expected. 

In figure [5l the temperature is plotted against p^. 
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Figure 5: This plots shows the behavior of the temperature as a function of . 

4.3 Entanglement Entropy 

The quantities H and [3 are: 

The difference in area and / are given as 

dy i y^- plyl 2 \ T* 



y* 



--dy 



(4.16) 



I = M djL_ (417) 

In the following graphs, A A is plotted in units of UqR^ against Iuq and p^. 

In figure [6], AA is plotted against for different values of here we observe that 
the difference in entropy is strictly bounded by zero from above. In figure [7|^b), AA is 
plotted against Iuq for values p^ less than 3. From this graph it is apparent that the 
difference in entropy approaches a constant negative value for large values Iuq. We also 
note that an increase in p^, translates to a decrease in AA. This feature is apparent 
in figure [6] and in figure [71 Furthermore, we also see that increasing y^, in figure [6l 
shifts AA downward. This allows us to also see that the difference in area is negative 
for large values of /; since small values of coincide with large values of / as shown 
in figure [7](a). From these observations, we conclude that the difference in entropy is 
always negative for p^ < 3. Thus there is no transition. 

For p^ > 3, we also do not observe any transition. This can be seen in figureOwhere 
the difference in entropy is strictly negative. In figure M^a) the difference in entropy is 
plotted against Iuq. Here we observe an increase in entropy for small p^ and then a 
monotonous decrease that is consistent with figure [61 
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Figure 6: The difference in entropy is plotted against at = Ip*, 2pl, Sp^,, 4p* (Blue, 
Green, Orange and Red). We observe that AA is always negative. 




012 4 6 8 10 1 2 3 4 5 



y* luo 

(a) (b) 

Figure 7: Figure [7|^a) shows the behavior of the order parameter, I, as a function of 
for = 0, 1, 2, 2.99 (Blue, Green, Orange, Brown). Figure[71|^b) shows the difference in 
entropy as a function of the order parameter. These plots show that there is no phase 
transition for p^ < 3 in the dyonic black hole asymptotic to ^^5*4. 

5 Global AdSp 
5.1 Global AdS^ 

Now we explore the entanglement entropy of black holes in AdS^ in global coordinates. 
This geometry does not satisfy the conditions of section [3] since there is no spacelike 
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(a) 



(b) 



Figure 8: Figure El^a) shows the behavior of / as a function of for = 
3.1,10,100,500,1000 (Brown, Blue, Green, Orange, Red). Figure EKb) Shows the dif- 
ference in entropy as a function of the order parameter. These plots show that there 
is no phase transition for > 3 in the dyonic black hole asymptotic to AdS4^. 

killing vector that commutes will with all killing vectors. We reformulate the problem 
explicitly in this case. The Schwarzschild black hole in global AdS^ is given by: 



ds^ = -fir)de + -rr^dr^ + r\d9^ + cos'fe)^^' + sin^f^) V) 
where the 3-sphere is written in Hopf coordinates and /(r) given by: 



f{r) 



1 - 



M 



L2' 



(5.1) 



(5.2) 



The conformal boundary is obtained, in these coordinates, by taking the r — > oo limit. 
The region A is given by the hypersurface on the boundary parametrized as: 

t = constant and < (p < 2ttI. (5.3) 

So, / measures the size of A. This quantity is bounded by 1 since the period of (j) is 2tt. 

5.1.1 Entanglement 

The surfaces in M that have the same boundary as A can be parametrized as: 

t = constant and r = r(0), (5.4) 
with boundary conditions as: 

r(0) = r(27r/) = r^, (5.5) 
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where is a cutoff parameter which will be taken to oo. The induced metric for this 
surface is given by: 



ds 



2 - + ^2 cos2(^))rf02 + r\de^ + sin2(^) V)- 



(5.6) 



Its area is then 



Ar = 2n I / rVcos2(e) + 62 sin(^) d9 d(j) where 6 = — ■)=f= (5.7) 
Jo 

2tvI rl 



= 2tx I I r^\/aP~+l^dad(f) where a = cos (6*). (5.8) 
Jo Jo 

The minimal surface is obtained by solving for the function r{(j)) that minimizes the 
action A^.. 

Ac is an action for a point particle with Lagrangian: 



1 = 271 f r^^/a? + 62 da. 
Jo 



From the Euler-Lagrange equation (H=Hamiltonian): 



(5.9) 



(5.10) 



We obtain the equation of motion 

d 







where c is some constant. Now we have 



dr' 



Wa^ + b^ 



Va2 + 62 



da 



' d 



da = —r^ I a—y/aP + l^da 
da 



Vr+62- / Va?+¥da 





(5.11) 



(5.12) 



(5.13) 



We can determine c by considering the point r^, where r' = 0. This corresponds to 
6 = 0. It is also important to note that must be bounded by tq. We thus have: 

c = -^. (5.14) 
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After integrating, we obtain the equation of motion of r(0) 

^3 



+ - ^2 In 



1 + Vl+F 



(5.15) 



We note that this is a transcendental equation for b. Thus we cannot evaluate the area 
explicitly. However we will be able to obtain some information when we substitute the 
equation of motion into the area formula. 

Before we proceed, we first evaluate the piece-wise smooth surface. This surface 
is the sum of 3 surfaces given by = 0, r = ro and = 27il at constant t. The line 
element for the (p =constant is, 

1 



fir] 



with area 



A = 271 



oo i-w/2 ^2 



>ro JO 

The r = ro piece has area 



sin(^) dOdr = 2it 







dr. 



l-n/2 

A = rj^(27r)2/ / sin 6 cos OdO = 27rVj^/. 
Jo 



The area of the piece- wise smooth surface is then 

^2 



Ah = 47r 



zdr + 27rVo/ 



(5.16) 



(5.17) 



(5.18) 



(5.19) 



5.1.2 Comparison 



Before we can compare the areas, we rewrite the area of the smooth surface as an 



integral over r. This is done by integrating out a; and then substituting for 
and integrated from r^, to r^o. One then obtains. 



Ar 



2-nl 

n I r^ 
'o 



Vl + 62 + 62 In 



2t\ 



Vl + 62 + 62 In 



l + v^TTF 



dr. 



We can also write / as a function of r*. This is done by integrating r'{(f)). 



2711 



2 — = 2 



r' Jr, hr^Ji^' 



(5.20) 
(5.21) 



(5.22) 
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Now we can evaluate the difference in area AA = — Act- It is given as: 



Air 



VlT¥ + In 



27r 



Air 



ry.2 



-dr - 2iirl 



dr 



dr 



(5.23) 



v^rT62 + ^2 In 



^dr. 



(5.24) 



Notice that we have spht the second integral in equation 15.131 to two integrals with 
ranges ro r* and r* ^ r^. Now we can use the equation of motion: 



i± = + - \? In 



1 + v^TTF 



(5.25) 



to obtain 
AA = 2ix 



'^oc ^2 



2¥ In 



2h 



* 



(ir — An 



-dr. 



Now we can take the r^, = ro limit and obtain: 
AA = Air 



^oo ^2 



dr < 



smce 



6 In I : I - 1 < 



for nonzero values of b. Thus we obtain 



AA(r, ^ ro) < 0. 



(5.26) 



(5.27) 



(5.28) 



(5.29) 



We can also evaluate the large r* limit by eliminating the In term in the difference 
equation: 



AA = Air 



Vl + &2 - 6 



rl + rl 
2r3 



dr — Air 



r, ^2 



ro 



--dr. 



(5.30) 
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As r^. — >■ Too, we have 



AA = -Air 



dr — An 



Thus we obtain 



2rb^ 



AA{r, Too) < 0. 



--dr. 



{5.31] 



(5.32) 



Now we see that the difference in entropy is negative for both large values of r^, and 
for r* = tq. We collect all this information in: 



AA(n = ro) < 0, and AA{r^ 



r„o) < 0. 



(5.33) 



From this argument, it is clear that if a transition exist, it must be at some intermediate 
value of r=K. Furthermore, if there is transition, there exist a r^: such that ^^^^ = 0. 
This may not be straightforward to do since A A is an integral equation over b{r,r^,)] 
which satisfies a transcendental equation. So one should integrate over b instead of r 
since the equation of motion can be read as: 



r{b,r^ 



VTT¥ - b^ In 



-1/3 



(5.34) 



One should also check that r has a one to one relation with b, which can be done 
by showing that the first derivative has no zeros. The derivative is then: 



dr 
db 



26(-l - VT+¥ + (1 + 62 + vO+F) ln( 



)) 



3vrTF(i + va+F)(va+F-62 in( 



))4/3 



(5.35) 



Now, db/dr is zero when the above expression diverges. However, this expression is 
always finite making it safe to integrate over b instead of r. The integration bounds 
for b will be from to oo. Notice that we have written A A in 3 different ways by using 
the equation of motion. In what follows, we use the one given in equation (15.301) . It is 



AA = in 



'/"oo ^2 



ViT¥-b 



2r3 



dr — Air 



ro 



-.dr. 



(5.36) 



It will also be useful to factor out the AdS radius L by making the coordinate change 
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X = r/ L and to introduce the parameter m = M/L"^. We then have: 
fir) 



M 
1- — + 



Xn 



J.2 IJl nf.2 

a/I + 4m - 1 



^ (x^ - xl){x^ + xl + l) = ^h{x, xo) 



x{b, y^) 

Ay4(x*,Xo) 



Vl+W - 6^ In 



1 + yiTF 



-1/3 



(5.37) 
(5.38) 

(5.39) 



X 



b^yh{x,Xo) 

^3 



2x3 



X 



xo 



\/h{x,XQ 



-.dx. 



Notice that we have introduced the polynomial function h[x,XQ) in (I5.37p . 
integral over b, the difference in area is given by: 



(5.40) 
As an 



AA 



oo ™4 „3 



xt9Hb)g'ib) 
by/h{x^g{b),xo 
x*{xl + xl)g'{b) 



[Vl + &2 - b]db 



X 



,0 v^/i(x,Xo) 



2b^yh{x^g{b),xo 



■.db where 



9{h) 



+ - 6^ In 



1 + v^TTF 



-1/3 



(5.41) 
(5.42) 



We can also write / as an integral over b: 



1 = ^ 



9'{h) 



-.db 



(5.43) 



TT 7o b^/h{x^g{b),xo) 

In figure [9] we have plotted the difference in area against / at various horizon posi- 
tions. The features observed here are similar to those in the planar systems. However, 
we see that the difference in area shifts downward as we increase the temperature in- 
stead of just scaling. Furthermore, we notice that there is a minimum r^, for which one 
can define the entanglement entropy. This minimum is given by the equation /(r*) = 1. 



5.2 Global ^^^4 
The relevant geometry is 

ds^ = -f{r)dt^ 



fir 



-dr^ + r^{de^ + sin" (9) d(j)^) 



(5.44) 
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(b) 



Figure 9: Figure [9t^a) shows the behavior of the order parameter, /, as a function of x* 
for xo = .5, 1, 1.5,2,2.5 (Brown, Red, Orange, Blue). Figure [9](b) Shows the difference 
in entropy as a function of the order parameter. These plots show no phase transition 
for the entanglement entropy in the Schwarzschild black hole in AdS^. 



where 



M 



(5.45) 



Introducing the x-coordinate as x = r/L, we find that the difference in area and 
the expression for / 

2L' I [x\q{b) - nb) - 2x1] - 2^^^ / .T^ (H .46) 

(5.47) 



AA 



b^/x^yh{x,Xo 

dx 



^ A'* by/x^Jh{x,XQ) 
The equation of motion and /i(x, xq) are given as 



37 X:jg 



^ where - ^"^^^ - ^'^(^' 



r'^f{r) x^f{x,Xo) 
h{x,xo) = x/(x, xo) = (x - xo)(x^ + XXo + Xq + 1). 

The position of the horizon is given by M as, 

2 M 
Xo + Xq - — = 0. 

The functions p{b) and q{b) are: 

q{b) = 

pip) = 



b^ + sm'^ {9)d9 



sm\9) 



y/b^ + sin^e) 



dO. 



(5.48) 
(5.49) 

(5.50) 

(5.51) 
(5.52) 
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By substituting the equation of motion in the difference of area integral, one can see 
that 

AA{x^ -^xo) <0 (5.53) 

since the integrand becomes negative for all b. Now we can plot the above function by 
integrating over b. In figure 15.21 we show the behavior of the difference in area with 
respect to /. Again we do not observe a transition. 




X* 1 

(a) (b) 

Figure 10: Figure [TO](a) shows the behavior of the order parameter, /, as a function of 
X* for xq = .5,1,1.5,2,2.5 (Brown, Red, Orange, Green, Blue). Figure [TOTb) Shows 
the difference in entropy as a function of the order parameter. These plots show no 
phase transition for the entanglement entropy in the Schwarzschild black hole in AdS^. 



6 Conclusions 

In lower dimensional systems the entanglement entropy has served as a way to detect 
exotic phases, see for example [29] . In the context of higher dimensional field theories, 
most interestingly, four-dimensional field theories, it was recently suggested that the en- 
tanglement entropy could be used as an order parameter for confinement / deconfinement 
[6]. Motivated largely by these developments we launched a systematic study of the 
entanglement entropy in what should be understood as field theories at finite temper- 
ature. The hope was that by studying the entanglement entropy we would be able to 
identify a transition, just as a study of the entanglement entropy in the supergravity 
dual of the confined phase suggested the possibility of a transition to the deconfined 
phase. The premier examples of supergravity backgrounds dual to field theories at 
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finite temperatures are the nonextremal Dp-brane solutions. We presented a detailed 
study in those geometries and generically detected no transition, except for the case of 
p = 6 which is of hmited interest from the field theory point of view. 

We then analyzed the entanglement entropy in global coordinates. The motiva- 
tion for the analysis comes from the fact that the Hawking-Page phase transition 
in asymptotically AdS spaces takes place only in global coordinates. This fact was 
given the interpretation, in the context of the AdS/CFT, of being dual to the confine- 
ment /deconfinement transition for a field theory on a sphere. We did not observe a 
transition for most of the supergravity backgrounds. Since, as shown in [16], the free 
energy undergoes a transition in this context and the entanglement entropy does not 
we can safely conclude that the entanglement entropy is not the free energy. This is 
an important observation since the free energy is considered the standard order pa- 
rameter for the confinement/deconfinement transition. The main reason arising from 
the natural expectation that in the confined phase the number of degrees of freedom 
comes from hadronic states and it is therefore of order for an SU{N) theory while 
in the deconfined phase one expects A^^ degrees of freedom. 

An interesting system where a Hawking-Page phase transition was established in 
[30] is the black hole with cascading asymptotics constructed in [31]. This system is the 
high temperature to which the Klebanov-Strassler background [32] flows. Preliminarily, 
with an analysis based on a few values of the temperature, we have been unable to 
observed a transition in the entanglement entropy for the background of [31]. We differ 
a more detailed analysis to a separate publication. 

The relationship between the entanglement entropy and the black hole entropy has 
been considered in various works in the literature, for example [34,35,37,36]. In [37] the 
one-loop corrections due to matter fields to the black hole entropy (computed through 
the free energy of the system) and to the entropy of entanglement were computed 
and found to agree for spin zero and spin one-half fields. A discussion more tuned to 
our main interest was presented in [34], where it was suggested that the microstates 
responsible for the black hole entropy are those due to the entanglement of the vacuum 
of the black hole. More bluntly, it was suggested in [34] that the leading contribution 
to the black hole entropy is due to entanglement. A similar proposal, albeit with a 
major modification about the position of the black hole, was argued in [35]. The idea 
that black hole entropy is due to the entanglement entropy, in fact finds a place in our 
analysis and was also explicitly realized already in [4] in the context of the holographic 
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dual of the A/" = 4 plasma. We believe this is, ultimately, the right track but, as we 
have shown in various examples, the story is more complicated due to the fact that 
the entanglement entropy does not behave as an entropy obtained directly from the 
free energy. It is worth stressing the significant role of the horizon, in fact, if one does 
not include the contribution to the piece-wise smooth configuration coming directly 
from the horizon, there is generically a transition. So, we conclude that the horizon is 
crucial in the thermodynamical competition. 

From a slightly more technical point of view, we can arrive at an important obser- 
vation regarding collapsing cycles. In all these geometries, as in the cases discussed 
in [6], there is a collapsing cycle coming from the temporal direction in Euclidean frame. 
However this cycle does not contribute in a significant way to the computations; in [6] 
the collapsing cycle was critical in obtaining transitions. An interesting alternative to 
the entanglement entropy has been recently suggested in [38]. The geometric entropy 
defined in this paper uses a double Wick rotation and therefore includes the full effect 
of the collapsing Euclidean cycle that was, in the original geometry, the temporal di- 
rection. We predict that this modification will lead to transitions in most of the cases 
presented in this paper. 

Finally, we end with a remark about the place of the entanglement entropy. We 
believe it is a thermodynamical potential rather than an entropy. Large parts of this 
intuition has been developed while collaborating with Carlos Nunez in similar issues. 
Namely, while studying the entanglement entropy, we observe that / and (the mini- 
mum of the smooth surface) behave like thermodynamically conjugate variables. This 
statement is motivated by the fact that when the function l{r^,) is double valued, the 
entropy exhibits branches as observed in the geometries in [6] and for p = 6. This 
would suggest that the entropy is a new thermodynamical potential. This is subject 
to further studies. Since we have observed, in a wide variety of models, no transition 
in the entanglement entropy we conclude that the entanglement entropy is, possibly,a 
thermodynamical potential that is different from the free energy. We are currently 
pursuing this line of investigation in collaboration with C. Nunez and expect to report 
on our findings soon. 
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